Necklaces with interacting beads: isoperimetric problems 
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Abstract. We discuss a pair of isoperimetric problems which at a glance seem 
to be unrelated. The first one is classical: one places N identical point charges 
at a closed curve T at the same arc-length distances and asks about the energy 
minimum, i.e. which shape does the loop take if left by itself. The second 
problem comes from quantum mechanics: we take a Schrodinger operator in 
L 2 (R d ), d = 2, 3, with N identical point interaction placed at a loop in the 
described way, and ask about the configuration which maximizes the ground 
state energy. We reduce both of them to geometric inequalities which involve 
chords of T; it will be shown that a sharp local extremum is in both cases 
reached by T in the form of a regular (planar) polygon and that such a T 
solves the two problems also globally. 



1. Introduction 

Isoperimetric problems are certainly one of the topics which keeps appearing in 
mathematical physics regularly. In the present paper we are going to discuss two 
new examples of this kind. It may seem that these problems differ mathematically 
and also have a different physical background; we will show nevertheless that they 
can be reduced to the same geometric question. A common feature is that they 
both concern extremal properties of interaction between N > 2 points placed at 
identical arc-length distances along a closed curve T of a fixed length L > 0. 

The first problem comes from classical electrostatics and concerns charged neck- 
laces. Let T : [0, L] — > M 3 be such a loop and suppose that N identical charges are 
placed at the points T(kL/N), k = 0, 1, . . . , N — 1. We ask about the shape which 
this constrained family of point sources will take in the absence of external forces, 
i.e. about minimum of the potential energy of the Coulombic repulsion. 

The other problem comes from quantum mechanics. Having again a loop T, now 
in two or three dimensions, we consider a class of singular Schrodinger operators in 
L 2 (M. d ), d = 2,3, which are given formally by the expression 

N-l , , 

(1-1) ff& = -A + fi£*(*- r (^ 

4=0 ^ ^ 
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We will recall below how they can be defined properly, see AGHH ; following the 
terminology of this monograph we can label the problem as a polymer loop. This 
time we are interested in the shape of T which maximizes the ground state energy 
of course, provided the discrete spectrum of H^ r is non-empty. 

Our goal in this paper is to show that the two problems reduce essentially to the 
same geometric question and that a sharp local extremum is in both cases reached 
by the shape with a maximum symmetry, in other words by a regular planar polygon 
with N vertices denoted in the following as Vn- Furthermore, we will show that 
the regular polygon represents also a global solution to the problem by reducing 
the task to a norm estimate of a particular operator on £ 2 (Z). 

The indicated quantum-mechanical isoperimetric problem was formulated first 
in the paper |Exlj , to which we refer for a deeper motivation, in the particular case 
when the curve T was an equilateral polygon. The problem was restated there in 
purely geometric terms and the existence of the maximizcr was proved locally. In 
|Ex2j a similar local result was derived for a continuous analogue of this problem in 
which the interaction was supported by the entire curve T and the maximizer is a 
circle. The technical improvement in the last named paper was that the geometric 
problem was viewed there from a general perspective in terms of inequalities for 
chords of the loop, more specifically, £ p norms related to the functions r(-+u)— T(-). 
In the subsequent work EHL a simple Fourier analysis was used to show that in 
the "continuous" quantum-mechanical problem the circle is a global maximizer. 

Our first aim here is to show first that the local proof of |Exl) can be extended 
to the more general class of curves. Then we will discuss how the argument of 
EHL can be modified to the present "discrete" situation. It appears that the task 
is more involved than in the "continuous" case, however, we will be able to reduce 
the question whether a regular polygon Vn represents a global extremum in both 
the isoperimetric problems described above to a well-defined operator problem. 



2. Point interactions on a loop 

To begin with, let us make more precise our requirements on the curve regu- 
larity. In what follows, we will suppose that for a fixed dimension d > 2 

(I) r : [0,L] — » M. d is a continuous, piecewise C 1 function such that T(0) = T(L) 
and |r(s)| = 1 holds for any s e [0,L] for which T(s) exists. 

The arc-length parametrization means in fact that we regard the curve as a map 
K — > R d (mod L). A shift in the argument is a trivial reparametrization, so without 
loss of generality we may assume that the point interactions are placed at 

(2.1) Vj :=T0£) , j = 0,l,...,iV-l; 

the indices can be again regarded as integers, yj = yj( mo dN)- A distinguished 
element of the described class is a regular polygon for which the points yj lie in a 

plane C K d (this is trivial if d = 2) at a circle of radius (2 sin . 

Let Yp := {yj : j = 0, ...,N — 1} be the interaction support. The object 
of our interest will the Hamiltonian — A Q .y r in L 2 (R d ) with N point interactions, 
all of the same coupling constant a 6 M. It is defined conventionally through 
boundary conditions which relate the generalized boundary values at each site yk, 
the coefficient at the singularity (logarithmic for d — 2, pole for d = 3) and the 
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next term in the expansion - see |AGHH] for a detailed discussion and recall that 
the construction of a point-interaction Hamiltonian does not work for d > 4. Recall 
also that a differs from the formal coupling constant a in ijl.lfl : it is sufficient to 
realize that the absence of a point interaction corresponds to a = oo. 

It is obvious that that spectral properties of — A Q .y r and — A Q .y r , correspond- 
ing to a pair of loops T and T' related mutually by Euclidean transformations of 
M. d are the same. This defines an equivalence relation on the set of all loops; with 
a terminological abuse we will speak about curves having in mind such equiva- 
lence classes. We will suppose that the operator — A Q ,y r has a non-empty discrete 
spectrum, 

(2.2) e x =ei(a,y r ) := infer (-A Q ,y r ) < , 

which is true for any a € M if d = 2, while in the case d = 3 it is a nontrivial 
assumption satisfied below a certain critical value of a - cf. jAGHHl Sec. II. 1]. 

One of the main results of this paper referring to the second of the problems 
mentioned in the introduction, the polymer loop, can be formulated as follows: 

Theorem 2.1. Assume (£) and if then ei(a, Yr) is for a fixed a and L > 
locally sharply maximized by a regular polygon, T = Vn ■ 

The proof will be done in several steps; in this section we will demonstrated that 
the task can be reduced to a purely geometric problem. With the usual notation, 
k = in with k > 0, we find the eigenvalues — n 2 from the spectral condition, 

det Q k = with (Q fe ) - := (a - £ fe )<^ - (1 - ^)g% 



where g& := Gk(yt — Vj), or equivalcntly 



( 2 - 3 ) Sij = \ ■ 



- Vj\) ■■■ d = 2 
rf = 3 



I *7r|j/j— Vj\ 

and the regularized Green's function at the interaction site is 

f -2^( ln f+7E) ... d=2 
I ... d = 3 



where 7e is the Euler number. The matrix Qi K has N eigenvalues counting multi- 
plicity which are decreasing in (— oo,0) as functions of k by jKLL IAGHH] . The 
spectral threshold e\(a, T) corresponds to the point n where the lowest of the indi- 
cated eigenvalues vanishes. Consequently, we have to check that 

(2.4) mma(Q ikl ) < mm a(Q lkl ) 

holds locally for T ^ Vn- Here and in the following the tilded quantities correspond 
always to regular polygon, T = Vn, in particular, —k\ = ei(a, Vn)- 

Next we use the fact that the lowest eigenvalue of Qik x corresponds to the 
eigenvector <j>i = A r_1/ ' 2 (1, . . . , 1), because by |AGHH] there is a one-to-one corre- 
spondence between an eigenfunction c = (ci, . . . , Cjy) of Qi K at the point, where the 
corresponding eigenvalue vanishes, and the corresponding eigenfunction of -A 0i y r 
given by c <-> J2j=i c jGi K {- ~ Vj), up to a normalization. Again by jACHHj . the 
principal eigenvalue of — A ai Y r is simple, so it has to be associated with a one- 
dimensional representation of the corresponding discrete symmetry group of Vn; it 
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follows that the coefficients are the same, c\ = ■ ■ ■ = cn- This yields the expression 



(2.5) 



min a{Q ikl ) 



^i) = «-e* i -^£<r- 



l<3 



On the other hand, for the left-hand side of l|2.4[) we have a variational estimate, 
min a(Q m ) < fa , Qmfa) = a - - A £ , 



N ■ 

i<j 

which shows that it is sufficient to check validity of the inequality 

(2.6) £ G «G/* - Vj) > £ Gin(Vi - Vj) 

i<j i<j 

for all k > and T ^ Vn- Let us introduce the symbol £ij for the chord length 
\Vi - yj\ and l VJ ;= \y l - jf,-|, and define the function F : (R + ) N( - N - 3 ^ 2 — > R by 

[N/2] 

f({^})--=J2 E [^(%)-gu4-)] ; 

m=2 |j— j\=m 

we want to show that F({£ij}) > except if {£ij} — {£ij}- The function Gi K (-) is 
by l|2.3|l convex for a fixed k > and d — 2, 3, thus by Jensen's inequality we have 



[N/2] 

F^}) > $>„ 

m—2 



G, 



*" i - -i 
|* — J \— m 

where v n is the number of the appropriate chords, 

j N ... m = l,...,[§(AT-l)] 

m = iN for N even 



±7V 

2 IV 



At the same time, Gi K {-) is monotonously decreasing in (0, oo), so the sought claim 
would follow if we demonstrate the inequality 



(2.7) 



1,-m+l > — ^i. 



I*— Jl 



and show that it is sharp for at least one value of m if V ^ Vn ', in this way we have 
accomplished the goal to reformulate the problem in geometric terms. 



3. Inequalities for chord sums 

It is useful to discuss the geometric problem without dimensional restrictions, 
i.e. for any d > 2. Denoting again for a given T by yj € M. d the points on the loop 
defined by i|2.1[l . we will study the following family of inequalities 

(3.1) Dl N (m) : £* =1 \y n+m y n \*> < , P > , 

(3.2) Dl%{m) : En =1 |j/n+m - S/nr > City , P > , 

for any m = 1, . . ., [|-/V], where [•] denotes as usual the entire part; they contain as 
a particular case the inequalities for equilateral polygons discussed in Ex2 . 
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The inequality (|2.7|) is nothing else than D\ N (m). In the next section we are 

going to demonstrate that it holds locally, i.e. in the vicinity of Vn, proving thus 
Theorem 12 .11 Furthermore, in Section [5] we will show that a stronger claim can be 
made, namely that the inequalities and 1)3. 2(1 hold generally as long as p < 2. 
Without imposing a restriction on TV this seems to be an optimal result, because 
the example of a rhomboid shows that D P L 4 (2) cannot be valid for p > 2. 

Notice that to establish the above described property one needs in fact to prove 
the inequality D\ N (m) only as the following simple result shows. 

Lemma 3.1. D P L N (m) implies D P L N (m) if p > p' > 0. Similarly, D P L N {m) 
implies Z)^^ v (m) for any p > 0. 



Proof. The first claim follows from convexity of x i— > x a in (0, oo) for a > 1, 
rl -PL p sin p — JL / ,\p/p' (^ ,\ P/P 



sin p £ ^ . . , 

N n=l \ n=l / 

it is then sufficient to take both sides to the power p' ' jp. On the other hand, 



N 2 N 1+p sin p $ 



iY 

E, _ |-p ^ > iv ' aiIr IV 

Wn+m Dn\ d_ N d. Tp-pxm 

holds by Schwarz inequality; it gives the second claim and completes the proof. □ 



4. Local extrema 

First we are going to prove Theorem 12. II As we have said one has to establish 
that the inequality l|2.7|l . or D\ N {m), holds locally. Since the argument repre- 
sents an extension of the proof of Theorem 4.1 in |Exlj to a more general class of 
constraints, we will skip details in repeating the latter. We are looking for a local 
maximum of the function of Nd variables, 

1 N 

fm ■ fm(Vl, ■ ■ ■ , Vn) = E ~ Vi+ri\ , 

1=1 

under the inequality-type constraints gi(yi, ■ ■ ■ , y n ) > 0, where 

9i(yi, ■■■>y n ) ■■= - \Vi~ Vi+i\, i=l,...,N; 

the true number of independent variables is (N — 2) (d — 1) — 1 because 2d — 1 
parameters are related to Euclidean transformations and can be fixed. 

Following the convention for inequality-type constraints we introduce slack vari- 
ables z r , r = 1, . . . , N, and Lagrange multipliers X r , r = 1, . . . , N, which determine 

JV 

(4.1) K m (yi, ...,vn,zi,.. .,z N ) := f m (yi, . . .,yjv)+2j Ar (drivi, ••■,Vn)-zl) ■ 

r=l 

The first thing to compute are the derivatives d Vj K m which are equal to 

1 (y 3 -y J+m + y.-y^X _ // f ., _ yj - Vj - X 

N [\vj-yj+m\ \yj-yj-m\j l ■ l 
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Without loss of generality we may consider a regular polygon in the plane spanned 
by the first two coordinate axes and to parametrize its vertices in the following way, 



y± m 



^( ± E cos ^( 2 -+ 1 )'E sin ^( 2 " + 1 



n=0 



N 



n=0 



N 



so that 



yj±m\ — 



E c ° s jj( 2n + 1 



fm-\ 

E sin ^( 2?i+1 ) 

Vn=0 



Hence the gradient components d yj K m will vanish for j = 1, 
choose all the Lagrange multipliers in Ij4.1|) equal to 



AT ' 
Af provided we 



(4.2) 



A = 



AfT, 



with 



E:=o 1 sin^(2n + l) 



A' 



sm 



2 7rm 



sin 2 f 



notice that this quantity is always nonzero. At the same time, one has to require 
vanishing of the derivatives 



d Zj K m — 2Xj. 



■3 ' 



3 = h ■ 



,N, 



which means that at the extremum all the slack variables vanish, Zj = 0. This is 
not surprising; one naturally expects critical points of the function f rn to be reached 
under given constraints with the neighbor distances maximal, i.e. for a polygon. 

From this point on the argument proceeds as in [Exlj . Evaluating the Hessian 
at the stationary point, one can reduce the question about its negative definiteness 
to verification of the inequalities 



7rm . 7ir 
(4.3) sm _ S m-> 



7r . nmr 
sm — sm 

N N 



2 < r < m < 



1 



-N 



or cquivalently, the inequalities U m -i (cos -j^) > \U,„-i (cos ^) | for Chebyshev 
polynomials of the second kind, which can be done directly. 

The obtained result provides also a local solution to our electrostatic problem. 

Theorem 4.1. Under the assumption (£) the Coulomb energy of a charged 
necklace is locally sharply minimized by a regular planar polygon, T = Vn ■ 

Proof. For a given nonzero charge q the potential energy equals 



9 2 E^- 



Vk\ 



? 2 E^E 

m—1 n—1 



Vn-\ 



Vn\ 



and since by Lemma |3. II the inequality D\ N {m) implies D L 1 N (m), the sum of all 
repulsion-energy terms is locally sharply minimized by Vn- D 



5. Global validity of mean-chord inequalities 

Let us look now what is needed to prove global validity of the inequalities (|3.1|) 
and l|3.2fl . in particular, to see whether a Fourier analysis in the spirit of EHL 
could help; in view of Lemma 13.11 it is enough to consider the case p — 2 only. 
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With the natural scaling properties in mind we may without loss of generality put 
L = 2-7T and to express the function T through its Fourier series, 

(5.i) r( s )= J2 c « ems 

0#n£Z 

with c„ G C d ; since T(s) G M. d the coefficients have to satisfy the condition 

The absence of the coefficient Cq means naturally no restriction; it can be always 
achieved by a choice of the coordinate system. 

It is convenient to impose slightly stronger regularity requirements on T as- 
suming that it is of the C 2 class; recall that validity of D 2 L N (m) can be extended 
from such a family of loops to those satisfying the hypothesis (I) by means of the 
Weierstrass theorem and continuity of the functionals involved. In such a case KF 
Sec. VIII. 1.2] the derivative of T is a sum of the uniformly convergent Fourier series 



(5.2) 



t(s) 



0#nGZ 



The assumed arc-length parametrization means |r(s)| = 1 giving thus the relation 
2n= \T{s)\ 2 ds = J2 E riltf -Cne^-^ds, 

where c[ = (ci t %, . . . ,ci y d) and dot marks the inner product in C d , or equivalently 
(5.3) 



E 



2 1 2 -i 

n c„ =1 



^-2-nimj/N _ ^ ^2-n 



imk/N 



1 e 



2-nin(k~j)/N 



Furthermore, using l|5.1|) we can rewrite the left-hand side of Z)| jv( to ) as 

N 

E E 

n=l o^j.fcez 

Next we change the order of summation and observe that X^n=i e 2mn ^ k "^/ N = N 
if j = k (mod N) and zero otherwise; this allows us to write the last expression as 

(5.4) 



4Ar E E IjX-i*i<* 



iez o =£ j, k e z 

j - k = lN 



j sm „ r 
AT 



A; 1 sin ■ 



If D 2 ^ ^(m) should be valid, this quantity must not exceed the right-hand side of 
(I3.1|) for p = 2; hence the sought inequality is equivalent to 



(5.5) 



d, (A (Ar,m) ® J)d) < 



7T sm 



AT 



TV sin jj 



where the vector d G £ 2 (Z) £g) C d has the components dj :— \j\cj and A( N,m ) is an 
operator on £ 2 (Z) determined by its matrix representation as 
(5.6) 



.4 



(N,m) 
jk 



N 



if ^ j,k eZ, j = k (mod AT) 
otherwise 
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It is obvious that A^ N,m ^ is bounded because its Hilbert-Schmidt norm is finite. 
Since d is a unit vector in view of 1)5.3)1 we arrive at the following conclusion: the 
inequality D\ N (m) with fixed N, m is valid provided the norm of the operator 
A( N ? m ) does not exceed the quantity at the right-hand side of l|5.5|l . 
This sufficient condition allows us to prove our second main result. 

Theorem 5.1. Within the class of curves specified by the assumption (£) the 
inequalities D 2 L N {m), and thus also $3.1]) and with p < 2 for fixed values of 

N = 2,3, . . . and m = 1 , . . . , , are valid. 

Proof. For a given j : =/= and d e £ 2 (Z) the relation (|5.6(l gives 

nmk 



_ i . irmj 
7 sin -—— 

N 



E 



o / fe e z 

k = j (mod N) 



k- 1 



sm ■ 



N 



The norm H^l^'"^^! is then easily estimated by means of Schwarz inequality, 



2 nm J 
N 



E 



o ^ k e z 

k = j'(mod N) 



i irmk 




k sm 


dk 


N 





N-l 

< e sin4 

n=0 



N 



S ... 



E 



1n+lN\ 



n + IN ^ 

I e Z 



where we have introduced 
S n := E 



f 



n + IN =f 
I 6 Z 



{n + lNf 



E 

/=i 



(IN-n) 2 {IN-N + n) 2 



However, the above series is easily evaluated to be 

2 

/ 7T 1 

Sn 



N sin s» 



and since ||<i|| 2 = X) n =o Siez l^n+HV | 2 > the sought inequality follows from 1)4.3)1 . □ 
This allows us to strengthen our claims concerning the two original problems. 

Corollary 5.2. Adopt the assumptions (£) and 1)2. 2f) : then ei(a, Yr) is for a 
fixed a and L > globally maximized by a regular polygon, T = 7-V- 

Corollary 5.3. Under the assumption ^ the Coulomb energy of a charged 
necklace is globally minimized by T = Vn- 

Remark 5.4. Notice that the inequality 1)4.3)1 used in the proof is sharp. That 
means that, in distinction to the "continuous" analogue of our problem, the ex- 
tremum cannot be reached in the class of C 2 smooth functions in which we have 
performed the described Fourier analysis. 
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6. Concluding remarks 

Let us first comment on relations to the "continuous" case treated in EHL 
where the global validity of the inequalities analogous to 13. 1|) and t\'d. 2f> was proved. 
Notice that formally that situation corresponds to N = oo. The counterpart of the 
operator A^ N ' m ^ is then a multiple of the unit operator and it is only necessary to 
employ the inequality lsin^p-1 < \j sin 2 ^], or slightly more generally 

(6.1) | sin j a; | < j sinx 

for any jgN and x £ (0, ^ir], which is checked easily by induction. In the present 
case with a finite N the operator has infinitely many side diagonals such a simple 
estimate based on (|6.1|) is too rough, because it yields an unbounded Toeplitz- 
type operator, and one has to do better using the matrix-element decay in 1(5.6(1 . 
Fortunately it can be done as the proof of Theorem 15 . 1 1 shows . 

Another aspect of the relation between the two cases is that in the continuous 
case the analogue of (|5.4|) follows from Parseval relation and the quantity is natu- 
rally invariant with respect to shifts in the arc-length parametrization. This is not 
true here; recall that the shift s — > s + sq is equivalent to the replacement of Cj by 
Cje is ° , which changes in general l|5.4(l due to the presence of the off-diagonal terms. 

Let us also comment briefly on various extensions of the present problem re- 
stricting ourselves to the "discrete" situation only. A natural question concerns 
the existence and properties of the extrema in situations when we have no built- 
in symmetry, either by assuming a nonconstant sequence of coupling parameters 
{<Xj} and/or taking their sites Sj at the loop with a non-equidistant distribution. 
In both cases the task becomes more difficult because we can no longer use the 
relation 1(2.5(1 which lead us to the geometric reformulation based on the inequality 
(12.6(1 . in particular, the solutions for the charged necklace and polymer loops will be 
now in general different. Another extension could concern point interaction family 
in M 3 placed on a closed surface. It is again not straightforward, however, because 
in distinction to a curve such a surface cannot be locally rectified and the answer 
will depend at the choice of the source sites at the surface. 
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